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O '' Abstract 



The effects of absence of inversion symmetry on superconducting states are investigated the- 
oretically. In particular we focus on the noncentrosymmetric compounds which have the cubic 
Q ■ symmetry O like Li2Pt3B. An appropriate and isotropic spin-orbital interaction is added in the 

o : 

J . Hamiltonian and it acts like a magnetic monopole in the momentum space. The consequent pairing 

^ \ wavefunction has an additional triplet component in the pseudospin space, and a Zeeman mag- 



netic field B can induce a collinear supercurrent J with a coefficient k{T). The effects of anisotropy 
embedded in the cubic symmetry and the nodal superconducting gap function on n(T) are also 
considered. From the macroscopic perspectives, the pair of mutually induced J and magnetization 
M can affect the distribution of magnetic field in such noncentrosymmetric superconductors, which 
is studied through solving the Maxwell equation in the Meissner geometry as well as the case of 
a single vortex line. In both cases, magnetic fields perpendicular to the external ones emerge as a 
signature of the broken symmetry. 
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I. INTRODUCTION 

The family of fermion superfluid, which includes the classes of conventional supercon- 
ductor, heliuni-3 superfluid, and cuprate superconductor, has been one of the most frontier 
subjects in condensed matter physics. According to the parity symmetry of their pairing 
wavef unction,- the above classes can be labeled as s-wave, p-wave, and c/-wave superfluid 
respectively and each has distinct thermodynamic and transport properties. In a system 
without inversion symmetry, this classification is however invalid, and the system is ex- 
pected to simultaneously possess the properties belonging to distinct classes. Symmetry of 
the pairing wavefunctions as well as the gap functions are the immediate question. Theoreti- 
cal studies based on the addition of a spin-orbital interaction in the Hamiltonian predict the 
Cooper pair to be a mixed state of singlet and triplet in pseudo-spin space,- which can lead 
to a non-vanishing spin susceptibility at zero temperature.->^i^ Besides, the nodal gap struc- 
ture has been investigated experimentally^'^'^i^ on two typical superconducting compounds, 
CePtaSi and Li2Pt3B,->i^ which have the point group symmetries of C^v and O, respectively. 

On the other hand, the spin-orbital interaction also provides a correlation between the 
electric and magnetic degrees of freedom within the Fermi sea, connecting with the magnetic 
properties of the superconducting state in a subtle way. For example, a net polarization of 
spins can be induced by a shift of momenta distribution or vice versa in the superconducting 
state.->ii>i^ In other words, the supercurrent and magnetization can be mutually induced. 
Therefore, the macroscopic distributions of current and magnetic field in the superconducting 
state can also be used to probe the effect of lacking inversion symmetry.— li^ii^ii^ One should 
note that the the form of spin-orbital interaction must vary with the background crystals 
of different point group symmetry. Consequently the magnetic properties pertaining to 
superconducting CePtsSi and Li2Pt3B are expected to be quite different. However, almost 
all the previous theoretical studies are based on the symmetry of C4i, which allows the 
Rashba form of spin-orbital interaction. 

In this paper we focus on the magnetic properties pertaining to the compounds with 
crystal symmetry of O, such as Li2Pt3B. The starting point is to write down an appropriate 
spin-orbital interaction which turns out to act like a magnetic monopole in momentum space 
in this case. For simplicity we first consider the case of isotropic Fermi surface and pairing 
gap. The supercurrent J is found to have a component parallel to the applied Zeeman 



magnetic field B and the proportional constant n is obtained by the linear response theory. 
For macroscopic studies, we employ the Maxwell equation to investigate the distribution of 
magnetic field in the Meissner geometry and the case of a single vortex line. Lastly, we also 
consider the effects of anisotropy embedded in the cubic symmetry which causes a power- 
law dependence of k,{T) for very low temperature due to the appearance of line nodes of 
superconducting gap functions. 

II. MICROSCOPIC DERIVATION OF SUPERCURRENT INDUCED BY A ZEE- 
MAN FIELD 

The goal of this section is to obtain an expression for the supercurrent induced by a 
Zeeman magnetic field in the bulk superconductor without inversion symmetry. We first 
consider the normal state. The lack of inversion symmetry is manifested itself by the spin- 
orbital interaction in the Hamiltonian H = J2p{Hp)ai30'pa(^pf3' ^^ which the operator is given 
by, 
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Hp = ^-EF-K-a, (1) 

and apa are the second-quantized operators for electron of momentum p and spin polar- 
ization a = {ti} is along the z axis in the laboratory frame. For convenience we write 
,^p = p'^/2m — Ef. Note that the spin-orbital interaction is characterized by the parity- 
breaking inner product consisting of a parity-odd hp = — /i_p and the spin a which is 
invariant under spatial inversion. It is convenient to work in the helicity basis (labelled by 
f|~J|) in which the operator hp ■ a is diagonal, that is, 

hp-a\pU)=±\pn)- (2) 

The eigenvalues of Hp are thus given by Cp = Cp =F hp for positive fl^ and negative J| helicity, 
respectively. Hence the degenerate spectrum is split into two branches ± in the presence 
of the spin-orbital interaction. The transformation between the helicity basis fl^JJ- and the 
laboratory frame basis fl are given by the unitary operator Up = exp(— |k ■ aOp) which 
rotate the z axis by an angle of 6p around the axis of k = z x /ip. More explicitly the matrix 
form of Up can be written down in terms of the coordinate of hp, namely. 



f/n 



I 



COS ■ 



j^p sin % 



-i(f)p 



sm 



cos -Y 



(3) 



Next we include the pairing between two electrons of opposite momenta on the same 
branch. In the helicity basis, a general mean-field description for the pairing potential H/^ 

can be written as 



Ha = J2 [^+(p)«-pfr«pt + ^-(p)«-pJ|«pJ| + h.c. 



(4) 



where the two gap functions A+ and A_, representing the pairing order parameter on the 
two branches, are not identical in general. However, the above pairing Hamiltonian can, 
by performing the transformation U, be restored to the case of singlet whenever A-|-(p) = 
e=F*0p|/\| Now the Nambu representation for the full Hamiltonian H in the helicity basis 
can be written as 



h=j: 



+ clLh a 



%t ^-pt 



( 



ip-K 



A_ 



p^ "-P^ 



1^ a; -^p + K 



A*_ -e 



p '^p 




(5) 



In what follows we employ the method of Matsubara Green's functions. ^^ It is useful to 
introduce the Nambu spinor representations \E'p and \&p for 1|-I| and t| basis, respectively 



^. 



^ n ^ 

"pH- 

Op^j, 

«-pi> 



/ 



' ^p 



\ 



^Pi 



'-pT 



V «-Pi / 



(6) 



The Matsubara Green's functions G in the fj basis are defined in a complex time-ordered 
manner as 



G„;3(p,r) = -(r.vl>p„(r)v[/t^^(0)) 



(7) 



and in the matrix form is 



G 




(8) 



where g and / are the matrix forms of the ordinary Green's functions and Gor'kov Green's 
functions. We note that the lower components have the properties f7(p,T) = — 5'*'^(— p, — t) 
and /(p, r) = /^(p, — t). The Fourier transformation of G is given by 



G„^(p,cu„) = \ rrfre'-"^G,^(p,r) . (9) 

p JO 

where 1//5 is the temperature and the frequency cij„=(2n + l)7r//3 is restricted due to the 
Fermi statistics. It is easier to first obtain the Green's function G by simply inverting the 
matrix {iuJn — H) in the helicity basis. The desired G can be obtained by performing the 
rotation in the pseudospin space using the foUowings, 



/(P) 



ujui 



p/f/-p 



(10) 

(11) 



where the transformation matrix f/_p for the opposite momentum is given by Up{—i<7y)e 



iaz4>p 



Using the property that ayff ay 
obtained as follows, 



-a, the expression for the Green's function G can be 
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f 



2 

1 r 



i9+ + 9-) + i9+-9-)hp-^ 

(/+ + /-) + (/+ - f-)K ■ ^ 



Uo-„ 



where the scalar functions g± and f± are given below 



(12) 



iuJn + e 



± 



9± 



h 



^2+g±2^|A 12 ' 



A±e 



±id 



ujl + e^^ + \^ 



(13) 



where we note that the previous condition for the pairing to recover the singlet is consistent 
with the condition for which the triplet component of Gor'kov Green's function vanishes in 
Eq. ^. 



In what follows we use the linear response theory to calculate the supercurrent J in- 
duced by an external Zeeman magnetic field b. We express the Fourier-transformed current 
operator in terms of the Nambu spinor representation as 

p 

where the momentum p± = p ± |, and the velocity operator Vp associated with momentum 

p is obtained by taking derivative of Hamiltonian with respect to p, which gives an identical 
result with the previous studies,— 

Vp = " ' ' . . • (15) 

^ + Vp/ip - 

The paramagnetic perturbation V resulting from the Zeeman magnetic field b{r) is V^ = 
—fi J drs{r) ■ b{r), here the positive /i is the magnetic moment and s{r) denotes the local 
spin density. V can be represented in Fourier space as V^ = — /iSq ■ 6_q in which the Fourier- 
transformed spin density is given by 

^"'q = E<-q/2S^P+q/2, (16) 

P 

where S is the spin operator for the Nambu spinor representation, given by 



(17) 





After some arrangement the current J can be written down explicitly as: 



rf^ 



J{q,iu;n) = - f rfre*""^(Tj;(r)sq(0)) ■ (-/i6_q) 

^ U 

= -^EEoTr' VpG(p_,iu;m)(S-b_q)G(p+,ia;n + i^i, 



(18) 



in which the symbol Tr' denotes taking trace over both the electron and hole sectors, and a 
factor of I is added to avoid double counting. 



III. MANIFESTATION OF ABSENCE OF INVERSION SYMMETRY 

In this section we are going to demonstrate the manifestation of absence of inversion sym- 
metry in this cubic superconductors from both microscopic and macroscopic perspectives. 
Starting from the general expression for supercurrent in Eq. (1181) . we investigate the static 
and homogeneous case, that is, the hmits q — > and Un ^ are taken. The resultant static 
current J is coUinear with the applied field, which can be written as 

J = -k{T)B (19) 

where the macroscopic magnetic field, or the magnetic induction is B = Bq. The appearance 
of this coefficient k, is an important signature of the lack of inversion symmetry. In Sec JIII Al 
K is studied explicitly for a given isotropic hp. 

Sec JIII Bl is devoted to the studies of macroscopic aspects, which deals with the interaction 
between the magnetic field and the nonvanishing pairing order parameter A. A crucial 
addition of Eq. flT9|) to the ordinary supercurrent and the corresponding magnetization are 
the key ingredients for understanding the new distribution of magnetic field. By the way, 
the expressions for J and M can be obtained by taking the derivatives of free energy as,— 



J -2^ 

6q 

M=-^, (20) 

where the gauge-invariant phase gradient q = hVip + —A and the free energy F contains 
an extra term of — ^/tq- B representing the absence of inversion symmetry.— More explicitly 
the above expressions can be written as 



_eJ = _(A+-V^)--B, (21) 






47rM = —{A+—Vcp). (22) 



A denotes the London penetration length. The length parameter S = -^k\ is introduced 



c 



for later convenience. Note that we have taken the electronic charge to be (— e) 



A. Microscopic aspects 

Now we consider the simplest case for which the spin orbital interaction is isotropic, and 
the gaps are identical for both branches and isotropic as well, that is, 

— * 

hp = ap, 
A+ = A_ = A . (23) 

The strength of spin-orbital interaction is characterized by the quantity a, which has the 
dimension of velocity and is weak in the sense that a/vp ^ 1. Here vp denotes the Fermi 
velocity. Starting from Eq. (TT5]) . with some arrangements, the static current can be obtained 
by taking the limit q -^ of the following expression, 



where the matrix elements of Q represent the factors for intra-branch contributions, fi = 
u = ±, and the inter-branch ones, fi = —v. Explicitly, Q in the matrix form can be written 
as 

Q= ' 1 ^ h (25 

where the three vectors are obtained after the operation of trace and solid-angle integration, 
namely. 



n 



/^IV|P(.-.B)(p..-)| = |PB. 

J Att m Sm 

I = j — Tr[aa((T ■ B)] = 2aB , 
J An 

f d^ 2 

t = J ^Tr[«a(p ■a)ia- B)(p ■ a)] = --aB . (26) 



Note that the trace here is only taken over a two-by-two helicity space in contrast to the 
previous operation in Eq. flTSl) . First we note that the static current in Eq. fl24l) is zero 
when the spin-orbital interaction is absent. For n, / and t in Eq. (l26l) . therefore, only the 
contributions up to first order a/vp are relevant. Since the summation over Un will give a 
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singular integrant concentrated at the Fermi level, it is eligible to substitute the quantities p 
and v{S,) with their values at Fermi level and then move them out of the integral. However, 
the contributions from n in the diagonal parts of Q must be taken care of because explicitly 
they are of zeroth order of ajvp. Hence the implicit contributions from modification of 
Fermi momentum and density of states due to the spin-orbital interaction have to be taken 
into account. Namely, the Fermi momentum for each branch, pp = pi?(l ± a/vp), and also 
the density of states at Fermi level, z/^ = mpp{l ± a/vp), should be used here. We also use 
the trick which enables performing the integration of energy first,— and after some algebra 
the coefficient k is obtained as 



^'"' -4ll^-^ErT^l-|i-^E ,.,.\.„.,., ":, J .(^r) 



The term in first bracket is actually the Yoshida function Y{A,T), which is a universal 
function characterizing the single-particle excitation across the gap A at temperature T. 
For an isotropic gap and at low temperature T/A ^ 1, the function Y is proportional to 
exp(— A/T). The second bracket is identical to the first one when the spin-orbital interaction 
is absent. We denote this term by a function y{A, a, T) to represent the excitations between 
two superconducting states separated by an energy of app. Thus we can rewrite Eq. ( I2711 
as K = |/iQ;i/(0)(|/ — Y). At zero temperature, the function y can be evaluated by replacing 
the summation over the Matsubara frequency by an integral, and for isotropic A and a this 
function is given by, 

1 



y(A, a, T = 0) = 1 - ^^7^^|f In (l + 26' + 2ev/rT^) , (28) 

where the number 9 = app/ A. For small 9, the function y ~ ^9"^. For large 9, it is 
approximately 1 — gr^- Both limits coincide with the previous predictions. ^'^^ 

In addition, the intra-branch and inter-branch contributions can be respectively recog- 
nized as the Pauli and Van Vleck ones in the previous studies.— Therefore the induced 
current J is absent as the Pauli and Van Vleck contributions cancel each other in the nor- 
mal state in which Y = y = 1. On the other hand, the existence of such current relies on the 
fact that the Pauli paramagnetic contribution in the superconducting state is significantly 
suppressed while the Van Vleck one is only reduced by a small portion as long as A ^ app. 



Consequently one can easily infer that the net supercurrent always flows in opposite to the 
Pauli paramagnetic current. 

B. Macroscopic aspects 

Here the effect of lacking inversion symmetry on a macroscopic length scale is studied 
through solving the static Maxwell equation, 

Vx B = 47rVx M+— (-e)J, (29) 

c 

Together with the current and magnetization given by Eq. (!2T]) and fl22|) . we are able to 
obtain an equation in terms of magnetic flled B only, namely. 



Vx VxB = -— B + 2— VxB, (30) 



in which the last curl term is generated from V x M as well as the collinear supercurrent 
induced by B. Hence, one can expect to observe a transverse component of the applied 
Zeeman fleld in such noncentrosymmetric superconductors. 

Eq. (I30l) can be applied for studying the penetration of magnetic fleld in Meissner geom- 
etry. Explicitly, we can consider a cubic superconductor occupying the space ior z > 0. It 
is more convenient to flrst consider a general fleld Bx{z)x + By{z)f containing both x and y 
components. Consequently, the equation with which the general fleld satisfles is, 

(31) 

where -B+ stands for the linear combination B^ + iBy. Deflning B^{z = 0+) = -Bin,+, the 
fleld just inside the superconductor, the general solution is then given by. 



d^ , , 


1 , . 


b d 


rf.^^-(^) - 


- A^^-(^' - 


-2'p..«-(^ 



^ a2+*A 



5+(^ > 0) = A„,+e V^ /, (32) 

which is identical to the previous results.— So one can expect a slight increase of penetra- 



tion depth by a factor of l/\/l — p- for such cubic superconductors in Meissner geometry. 
Besides, we note that the additional oscillation is a consequence of a parallel component of 
J to B. 
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The unknown -Bin,+ in Eq. (l32l) can be determined from the boundary condition Bext = 
Bin — AttNLIz = 0+), which requires the knowledge of magnetization M, or equivalently the 
gauge-invariant q. In fact q can be obtained from integration of the relation {—i)B^{z) = 
dA^{z)/dz with given boundary condition at infinity. We can assume the homogeneity for 
phase (f throughout the superconductor, which is indicative of vanishing A at 2; = oo to 
ensure zero current there. Consequently the boundary condition at the surface can be shown 
to, up to first order of 6/X, have the following form, 

iASext,+ = An,+ • (33) 

A similar relation for B_ can be obtained from the above by taking complex conjugates on 
both sides. 

If the external field is Bgxt = Bx, then A^{0^) = as a result of Eq. (1331) . Consequently, 
^x(O^) = 0, which demonstrates that the parallel field is continuous across the surface. In 
fact, the magnetic field B^ and By inside the superconductor can be obtained by taking 
the real and imaginary parts of Eq. ( 132|) . respectively. Up to first order of S/X, the two 
components can be written as. 



B^ = B 



By = B 



6z 
cos ^, + 


6 
A 


sin 


6z 
A2 


6 
A 


Sz 
cos^, 


— 


sin 


6z 
A2 



e-^ , (34) 

e-i . (35) 



Along the direction of applied field, the field B^ penetrates into the superconductor with 
an additionally slow oscillation of period about A/|. On the other hand, My(0"'") is finite 
due to the existence of finite fiow velocity proportional to Ay at the interface. Hence a 
discontinuity for field By is generated at the interface. 



-'-'in,y -Dext,y " 



(36) 



Bcxt,x A 

which is different from the previous prediction for inversion-broken superconductor of C^^ 

symmetry where the discontinuity happens to the parallel field across the interface.— The 

functional form of By in Eq. fl35l) indicates that it has the largest magnitude jB at the 

surface, changes sign at z = X and then decays to zero while slowly oscillating. Furthermore, 
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the flux associated with the perpendicular By is zero. This is consistent with the conclusion 
drawn from Eq. fl33|) that Ax = at the interface since both A^{z = oo) and A^ are zero. 

Eq. (1501) can also be applied for studying a single vortex line as a macroscopic signature 
of lacking inversion symmetry. We consider the conventional case in which the vortex line 
is along the z axis, and the cylindrical coordinates are adopted here. The components of 
magnetic field are assumed to be B^{r) and B.^{r) along the directions of and z, respectively. 
The z and (j) components of Eq. fl30|) are given by, 



1 d , d , 1 

r dr dr A^ 

d .1 d . 1 

dr r dr A^ 



B.{r) 
BJr] 



-k-- — {rB^) 
Xr dr 

"-XTr^^ ' 



(37) 
(38) 



where we denote the dimensionless number k 
assume the following perturbation solutions. 



26/X for convenience. We can therefore 



B^ = Bf^ + R''Bf^ + ..., (39) 

B^ = RBf + K'Bf + .... (40) 

The zeroth order solution Bf^ of Eq. fl37|) is just the conventional single vortex line solution, 
given by -^^Ko^r/X) where $ is a quantum of flux ^ and Kq is the modified Bessel function 
of zeroth order. As can be seen in Eq. ( !38l) . now the transverse field B^ emerges as a result 
of the nonzero source proportional to Ki from the identity Kq = —Ki. Up to the first order 
of K, the transverse field can be written down in terms of the Green's function gi associated 
with Eq. fl38|) .^^ Defining the dimensionless variable x as r/A, it can be expressed as. 



Bi\x)l 



$ 



27rA2 



x dx gi{x , x' ) Ki{x') 



Ki{x) x'dx'Ii{x')Ki{x') + Ii{x) x'dx'[Ki{x')] 



„'M2 



(41) 



in which Ki and Ji are the modified Bessel functions of first order. The asymptotic behaviors 
of the transverse field distribution are. 



Bi}\x) 



flni 

2 X 



TTX —X 




(42) 
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Hence the transverse field increases from zero at tlie origin, readies its maximum at a distance 
of order A from tlie center and is followed by an exponential decay. The above extra magnetic 
fields non-collinear with the externally applied one can in principle be detected by observing 
the extra precession of polarized muons when their polarization is parallel to the external 
applied field. | for Li2Pt3B is of order 10"'^ using the spin-orbital splitting estimated by Lee 
and Pickett.^ 

IV. ANISOTROPIC FERMI SURFACE AND LINE NODES OF GAP 

In previous sections we demonstrate an induced supercurrent parallel to the external 
Zeeman field as a signature of lacking inversion symmetry in cubic superconductors. Actually 
the spin-orbital interaction appropriate for the point group O respects all but the elements 
connected to inversion in Oh- The odd-parity basis functions^ {pS^ + P^Y + P2Z; n = 1, 3, 5}, 
in which the cubic symmetry is embedded, belonging to Aiu representation within Oh still 
can be used to construct the vectors h-p. Similarly, the general gap function respecting the 
cubic symmetry, 

A(p) = (Ao(p) + dp-(T)(za,), (43) 

can have the component Ao(p) constructed from the even-parity basis functions belonging 
to Aig representation while the vector function dp, having identical symmetry properties of 
/ip, can be constructed from Aiu representation. Here both components can be nonzero since 
parity is no longer respected. An important feature of the gap function given in Eq. (H3!) is 
the possible appearance of zeros when the order parameters Ao(p) and dp can simultaneously 
be real after appropriate gauge transformations, which is true if the time-reversal symmetry 
is respected in the system. Consequently it is possible to realize the zeros of gap function 
when I dp I exceeds |Ao(p)| for some points on the Fermi surface. Gapless excitation can 
therefore exist in such superconductors by showing, for example, a power law temperature 
dependence of penetration depth.— 

In fact the nodal structure of gap function in the compound Li2Pt3B was shown to be line 
nodes through the observation of linear temperature dependence in penetration depth for 
very low temperature.-*^ Here we shall investigate the effects of anisotropy and line nodes on 
the coefficient k(T) near zero temperature. Since we are only interested in the regime of very 
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weak Zeeman field, the anisotropy, which could result in some nonlinear field dependence 
for stronger field regime, has little qualitative effect here. Hence only the line nodes of 
gap function are relevant to the low temperature behavior of k{T). The intra-branch, or 
the Pauli, contributions of k(T) in Eq. flTTI) for the isotropic case can be generalized here 
by directly replacing the gaps A± with |Ao| ± |d|. As for the inter-branch contribution, 
the gaps on both branches are in fact much smaller than the separation app in Li2Pt3B, 
which suggests little relevance of actual gap function this contribution. Moreover, even the 
appearance of zeros associated with the spin-orbital interaction a{Q) are also irrelevant to 
the T-dependence of this contribution at very low temperature as long as the zeros associated 
with a are not identical to those associated with the pairing gap. 

We thus define a dimensionless quantity 7(T) = ""^ '^^.^ ' to present the temperature 
dependence due to the line nodes at temperature close to zero. Furthermore, from previous 
arguments, only the intra-branch contributions associated with the pairing gap of |Ao| — |d| 
is significant here. In addition, we are only interested in the effects due to the line nodes and 
take these parameters a{Q) and vf{^) to be isotropic, which makes the evaluation easier 
and more accessible. Next the summation over Un in the function Y in Eq. ( 1271) can be 
transformed into an integral of energy, which makes 7(T) into the following form, 

„r)^/^(co.«)W_(n)f .^(Jj^ + ^^^ji^) . (44, 

where the Zeeman field is assumed to be along the z axis, and A_{Q) denotes the gap 



functions of the direction fl on the Fermi surface, and E = y^^ _|_ A^(ii). We note that 
the above integral vanishes for T is exactly zero since E is always positive for all ^. Hence 
the contribution for T slightly larger than zero comes from integration around the solid 
angles Q associated with the zeros of gap. By the cubic symmetry, we can infer that 
there are six sets of line nodes on the Fermi surface, which as a whole remain invariant 
under any cubic rotation. Hence the contributions from the six sets can be divided into 
7(T) = 27||(T) + 47_l(T), where 7|| denotes that from the line nodes that are symmetrically 
distributed around the z axis, while 7_l denotes remaining sets that are around the x or y 
axis. 

For a given set of line nodes, the gap function can be expanded around these zeros in the 
following manner. 
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A4^)=A'ie,){9-e,), (45) 

which means the sohd angel ^2 = {6, 0) associated with the zeros can be parameterized 
as {9c{(/)),(/)) along the azimuthal direction. A'(^c) denotes the slope of gap function along 
the direction of 6. Finally the linear temperature dependence can be extracted from 7(T), 
and the following expression can be obtained if one extends the upper limit of 6 integral to 
infinity. 



711 (T) 
7±(T) 



2 - • ' ^ 




. X dxdy —^—- + TT^ — ,2/ /ox ' (46) 

Jo \r^ e"^ + 1 2r"'cosh (r/2)/ 

in which x = 9A'{6c)/T, y = ^/T and r = \^x'^ + y'^. In terms of polar coordinate the 
integral is to give 7rln2. 

V. CONCLUSIONS 

In this work we demonstrate an induced supercurrent parallel to the external Zeeman 
magnetic field utilizing the Green's function method. Besides, the induced supercurrent and 
the consequent magnetization modify the distribution of magnetic fields in the Meissner 
geometry as well as in the vortex line. Transverse magnetic fields are generated as a sign 
of breaking inversion symmetry in superconductors of point group symmetry O such as 

Li2Pt3B. 
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